Improved Indicators for Subnational Unmet Antiretroviral Therapy Need in the Health System: Updates to the Naomi Model in 2023.

BACKGROUND: Quantifying subnational need for antiretroviral therapy (ART) for HIV is challenging because people living with HIV (PLHIV) access health facilities in areas that may differ from their residence. We defined and demonstrated new indicators for PLHIV treatment needed to guide health system target setting and resource allocation. SETTING: Botswana. METHODS: We extended Naomi, a Bayesian small-area model for estimating district-level HIV indicators from national household survey and HIV service delivery data. We used model outputs for ART seeking probabilities in neighboring districts to define the "PLHIV (attending)" indicator representing the estimated number of PLHIV who would seek treatment at health facilities in a district, and "Untreated PLHIV attending" representing gaps in ART service provision. Botswana 2021 district HIV estimates were used to demonstrate new outputs and assess the sensitivity to uncertainty in district population sizes. RESULTS: Across districts of Botswana, estimated adult ART coverage in December 2021 ranged 90%-96%. In the capital city Gaborone, there were 50,400 resident PLHIV and 64,200 receiving ART, of whom 24% (95% CI: 20 to 32) were estimated to reside in neighboring districts. Applying ART attendance probabilities gave a "PLHIV attending" denominator of 68,300 and unmet treatment need of 4100 adults (95% CI: 3000 to 5500) for Gaborone health facilities. The facility-based "PLHIV attending" denominator was less-sensitive to fluctuations in district population size assumptions. CONCLUSIONS: New indicators provided more consistent targets for HIV service provision, but are limited by ART data quality. This challenge will increase as treatment coverage reaches high levels and treatment gaps are smaller.

Spatio-temporal estimates of HIV risk group proportions for adolescent girls and young women across 13 priority countries in sub-Saharan Africa.

The Global AIDS Strategy 2021-2026 identifies adolescent girls and young women (AGYW) as a priority population for HIV prevention, and recommends differentiating intervention portfolios geographically based on local HIV incidence and individual risk behaviours. We estimated prevalence of HIV risk behaviours and associated HIV incidence at health district level among AGYW living in 13 countries in sub-Saharan Africa. We analysed 46 geospatially-referenced national household surveys conducted between 1999-2018 across 13 high HIV burden countries in sub-Saharan Africa. Female survey respondents aged 15-29 years were classified into four risk groups (not sexually active, cohabiting, non-regular or multiple partner[s] and female sex workers [FSW]) based on reported sexual behaviour. We used a Bayesian spatio-temporal multinomial regression model to estimate the proportion of AGYW in each risk group stratified by district, year, and five-year age group. Using subnational estimates of HIV prevalence and incidence produced by countries with support from UNAIDS, we estimated new HIV infections in each risk group by district and age group. We then assessed the efficiency of prioritising interventions according to risk group. Data consisted of 274,970 female survey respondents aged 15-29. Among women aged 20-29, cohabiting (63.1%) was more common in eastern Africa than non-regular or multiple partner(s) (21.3%), while in southern countries non-regular or multiple partner(s) (58.9%) were more common than cohabiting (23.4%). Risk group proportions varied substantially across age groups (65.9% of total variation explained), countries (20.9%), and between districts within each country (11.3%), but changed little over time (0.9%). Prioritisation based on behavioural risk, in combination with location- and age-based prioritisation, reduced the proportion of population required to be reached in order to find half of all expected new infections from 19.4% to 10.6%. FSW were 1.3% of the population but 10.6% of all expected new infections. Our risk group estimates provide data for HIV programmes to set targets and implement differentiated prevention strategies outlined in the Global AIDS Strategy. Successfully implementing this approach would result in more efficiently reaching substantially more of those at risk for infections.

Subnational HIV incidence trends in Malawi: large, heterogeneous declines across space.

The rate of new HIV infections globally has decreased substantially from its peak in the late 1990s, but the epidemic persists and remains highest in many countries in eastern and southern Africa. Previous research hypothesised that, as the epidemic recedes, it will become increasingly concentrated among sub-populations and geographic areas where transmission is the highest and that are least effectively reached by treatment and prevention services. However, empirical data on subnational HIV incidence trends is sparse, and the local transmission rates in the context of effective treatment scale-up are unknown. In this work, we developed a novel Bayesian spatio-temporal epidemic model to estimate adult HIV prevalence, incidence and treatment coverage at the district level in Malawi from 2010 through the end of 2021. We found that HIV incidence decreased in every district of Malawi between 2010 and 2021 but the rate of decline varied by area. National-level treatment coverage more than tripled between 2010 and 2021 and more than doubled in every district. Large increases in treatment coverage were associated with declines in HIV transmission, with 12 districts having incidence-prevalence ratios of 0.03 or less (a previously suggested threshold for epidemic control). Across districts, incidence varied more than HIV prevalence and ART coverage, suggesting that the epidemic is becoming increasingly spatially concentrated. Our results highlight the success of the Malawi HIV treatment programme over the past decade, with large improvements in treatment coverage leading to commensurate declines in incidence. More broadly, we demonstrate the utility of spatially resolved HIV modelling in generalized epidemic settings. By estimating temporal changes in key epidemic indicators at a relatively fine spatial resolution, we were able to directly assess, for the first time, whether the ART scaleup in Malawi resulted in spatial gaps or hotspots. Regular use of this type of analysis will allow HIV program managers to monitor the equity of their treatment and prevention programmes and their subnational progress towards epidemic control.

Naomi: a new modelling tool for estimating HIV epidemic indicators at the district level in sub-Saharan Africa.

INTRODUCTION: HIV planning requires granular estimates for the number of people living with HIV (PLHIV), antiretroviral treatment (ART) coverage and unmet need, and new HIV infections by district, or equivalent subnational administrative level. We developed a Bayesian small-area estimation model, called Naomi, to estimate these quantities stratified by subnational administrative units, sex, and five-year age groups. METHODS: Small-area regressions for HIV prevalence, ART coverage and HIV incidence were jointly calibrated using subnational household survey data on all three indicators, routine antenatal service delivery data on HIV prevalence and ART coverage among pregnant women, and service delivery data on the number of PLHIV receiving ART. Incidence was modelled by district-level HIV prevalence and ART coverage. Model outputs of counts and rates for each indicator were aggregated to multiple geographic and demographic stratifications of interest. The model was estimated in an empirical Bayes framework, furnishing probabilistic uncertainty ranges for all output indicators. Example results were presented using data from Malawi during 2016-2018. RESULTS: Adult HIV prevalence in September 2018 ranged from 3.2% to 17.1% across Malawi's districts and was higher in southern districts and in metropolitan areas. ART coverage was more homogenous, ranging from 75% to 82%. The largest number of PLHIV was among ages 35 to 39 for both women and men, while the most untreated PLHIV were among ages 25 to 29 for women and 30 to 34 for men. Relative uncertainty was larger for the untreated PLHIV than the number on ART or total PLHIV. Among clients receiving ART at facilities in Lilongwe city, an estimated 71% (95% CI, 61% to 79%) resided in Lilongwe city, 20% (14% to 27%) in Lilongwe district outside the metropolis, and 9% (6% to 12%) in neighbouring Dowa district. Thirty-eight percent (26% to 50%) of Lilongwe rural residents and 39% (27% to 50%) of Dowa residents received treatment at facilities in Lilongwe city. CONCLUSIONS: The Naomi model synthesizes multiple subnational data sources to furnish estimates of key indicators for HIV programme planning, resource allocation, and target setting. Further model development to meet evolving HIV policy priorities and programme need should be accompanied by continued strengthening and understanding of routine health system data.

Evaluating distributional regression strategies for modelling self-reported sexual age-mixing.

The age dynamics of sexual partnership formation determine patterns of sexually transmitted disease transmission and have long been a focus of researchers studying human immunodeficiency virus. Data on self-reported sexual partner age distributions are available from a variety of sources. We sought to explore statistical models that accurately predict the distribution of sexual partner ages over age and sex. We identified which probability distributions and outcome specifications best captured variation in partner age and quantified the benefits of modelling these data using distributional regression. We found that distributional regression with a sinh-arcsinh distribution replicated observed partner age distributions most accurately across three geographically diverse data sets. This framework can be extended with well-known hierarchical modelling tools and can help improve estimates of sexual age-mixing dynamics.

Global and National Burden of Diseases and Injuries Among Children and Adolescents Between 1990 and 2013: Findings From the Global Burden of Disease 2013 Study.

IMPORTANCE: The literature focuses on mortality among children younger than 5 years. Comparable information on nonfatal health outcomes among these children and the fatal and nonfatal burden of diseases and injuries among older children and adolescents is scarce. OBJECTIVE: To determine levels and trends in the fatal and nonfatal burden of diseases and injuries among younger children (aged <5 years), older children (aged 5-9 years), and adolescents (aged 10-19 years) between 1990 and 2013 in 188 countries from the Global Burden of Disease (GBD) 2013 study. EVIDENCE REVIEW: Data from vital registration, verbal autopsy studies, maternal and child death surveillance, and other sources covering 14,244 site-years (ie, years of cause of death data by geography) from 1980 through 2013 were used to estimate cause-specific mortality. Data from 35,620 epidemiological sources were used to estimate the prevalence of the diseases and sequelae in the GBD 2013 study. Cause-specific mortality for most causes was estimated using the Cause of Death Ensemble Model strategy. For some infectious diseases (eg, HIV infection/AIDS, measles, hepatitis B) where the disease process is complex or the cause of death data were insufficient or unavailable, we used natural history models. For most nonfatal health outcomes, DisMod-MR 2.0, a Bayesian metaregression tool, was used to meta-analyze the epidemiological data to generate prevalence estimates. FINDINGS: Of the 7.7 (95% uncertainty interval [UI], 7.4-8.1) million deaths among children and adolescents globally in 2013, 6.28 million occurred among younger children, 0.48 million among older children, and 0.97 million among adolescents. In 2013, the leading causes of death were lower respiratory tract infections among younger children (905.059 deaths; 95% UI, 810,304-998,125), diarrheal diseases among older children (38,325 deaths; 95% UI, 30,365-47,678), and road injuries among adolescents (115,186 deaths; 95% UI, 105,185-124,870). Iron deficiency anemia was the leading cause of years lived with disability among children and adolescents, affecting 619 (95% UI, 618-621) million in 2013. Large between-country variations exist in mortality from leading causes among children and adolescents. Countries with rapid declines in all-cause mortality between 1990 and 2013 also experienced large declines in most leading causes of death, whereas countries with the slowest declines had stagnant or increasing trends in the leading causes of death. In 2013, Nigeria had a 12% global share of deaths from lower respiratory tract infections and a 38% global share of deaths from malaria. India had 33% of the world's deaths from neonatal encephalopathy. Half of the world's diarrheal deaths among children and adolescents occurred in just 5 countries: India, Democratic Republic of the Congo, Pakistan, Nigeria, and Ethiopia. CONCLUSIONS AND RELEVANCE: Understanding the levels and trends of the leading causes of death and disability among children and adolescents is critical to guide investment and inform policies. Monitoring these trends over time is also key to understanding where interventions are having an impact. Proven interventions exist to prevent or treat the leading causes of unnecessary death and disability among children and adolescents. The findings presented here show that these are underused and give guidance to policy makers in countries where more attention is needed.

Under-5 mortality in 2851 Chinese counties, 1996-2012: a subnational assessment of achieving MDG 4 goals in China.

BACKGROUND: In the past two decades, the under-5 mortality rate in China has fallen substantially, but progress with regards to the Millennium Development Goal (MDG) 4 at the subnational level has not been quantified. We aimed to estimate under-5 mortality rates in mainland China for the years 1970 to 2012. METHODS: We estimated the under-5 mortality rate for 31 provinces in mainland China between 1970 and 2013 with data from censuses, surveys, surveillance sites, and disease surveillance points. We estimated under-5 mortality rates for 2851 counties in China from 1996 to 2012 with the reported child mortality numbers from the Annual Report System on Maternal and Child Health. We used a small area mortality estimation model, spatiotemporal smoothing, and Gaussian process regression to synthesise data and generate consistent provincial and county-level estimates. We compared progress at the county level with what was expected on the basis of income and educational attainment using an econometric model. We computed Gini coefficients to study the inequality of under-5 mortality rates across counties. FINDINGS: In 2012, the lowest provincial level under-5 mortality rate in China was about five per 1000 livebirths, lower than in Canada, New Zealand, and the USA. The highest provincial level under-5 mortality rate in China was higher than that of Bangladesh. 29 provinces achieved a decrease in under-5 mortality rates twice as fast as the MDG 4 target rate; only two provinces will not achieve MDG 4 by 2015. Although some counties in China have under-5 mortality rates similar to those in the most developed nations in 2012, some have similar rates to those recorded in Burkina Faso and Cameroon. Despite wide differences, the inter-county Gini coefficient has been decreasing. Improvement in maternal education and the economic boom have contributed to the fall in child mortality; more than 60% of the counties in China had rates of decline in under-5 mortality rates significantly faster than expected. Fast reduction in under-5 mortality rates have been recorded not only in the Han population, the dominant ethnic majority in China, but also in the minority populations. All top ten minority groups in terms of population sizes have experienced annual reductions in under-5 mortality rates faster than the MDG 4 target at 4.4%. INTERPRETATION: The reduction of under-5 mortality rates in China at the country, provincial, and county level is an extraordinary success story. Reductions of under-5 mortality rates faster than 8.8% (twice MDG 4 pace) are possible. Extremely rapid declines seem to be related to public policy in addition to socioeconomic progress. Lessons from successful counties should prove valuable for China to intensify efforts for those with unacceptably high under-5 mortality rates. FUNDING: National "Twelfth Five-Year" Plan for Science and Technology Support, National Health and Family Planning Commission of The People's Republic of China, Program for Changjiang Scholars and Innovative Research Team in University, the National Institute on Aging, and the Bill & Melinda Gates Foundation.

Global, regional, and national incidence and mortality for HIV, tuberculosis, and malaria during 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013.

BACKGROUND: The Millennium Declaration in 2000 brought special global attention to HIV, tuberculosis, and malaria through the formulation of Millennium Development Goal (MDG) 6. The Global Burden of Disease 2013 study provides a consistent and comprehensive approach to disease estimation for between 1990 and 2013, and an opportunity to assess whether accelerated progress has occured since the Millennium Declaration. METHODS: To estimate incidence and mortality for HIV, we used the UNAIDS Spectrum model appropriately modified based on a systematic review of available studies of mortality with and without antiretroviral therapy (ART). For concentrated epidemics, we calibrated Spectrum models to fit vital registration data corrected for misclassification of HIV deaths. In generalised epidemics, we minimised a loss function to select epidemic curves most consistent with prevalence data and demographic data for all-cause mortality. We analysed counterfactual scenarios for HIV to assess years of life saved through prevention of mother-to-child transmission (PMTCT) and ART. For tuberculosis, we analysed vital registration and verbal autopsy data to estimate mortality using cause of death ensemble modelling. We analysed data for corrected case-notifications, expert opinions on the case-detection rate, prevalence surveys, and estimated cause-specific mortality using Bayesian meta-regression to generate consistent trends in all parameters. We analysed malaria mortality and incidence using an updated cause of death database, a systematic analysis of verbal autopsy validation studies for malaria, and recent studies (2010-13) of incidence, drug resistance, and coverage of insecticide-treated bednets. FINDINGS: Globally in 2013, there were 1·8 million new HIV infections (95% uncertainty interval 1·7 million to 2·1 million), 29·2 million prevalent HIV cases (28·1 to 31·7), and 1·3 million HIV deaths (1·3 to 1·5). At the peak of the epidemic in 2005, HIV caused 1·7 million deaths (1·6 million to 1·9 million). Concentrated epidemics in Latin America and eastern Europe are substantially smaller than previously estimated. Through interventions including PMTCT and ART, 19·1 million life-years (16·6 million to 21·5 million) have been saved, 70·3% (65·4 to 76·1) in developing countries. From 2000 to 2011, the ratio of development assistance for health for HIV to years of life saved through intervention was US$4498 in developing countries. Including in HIV-positive individuals, all-form tuberculosis incidence was 7·5 million (7·4 million to 7·7 million), prevalence was 11·9 million (11·6 million to 12·2 million), and number of deaths was 1·4 million (1·3 million to 1·5 million) in 2013. In the same year and in only individuals who were HIV-negative, all-form tuberculosis incidence was 7·1 million (6·9 million to 7·3 million), prevalence was 11·2 million (10·8 million to 11·6 million), and number of deaths was 1·3 million (1·2 million to 1·4 million). Annualised rates of change (ARC) for incidence, prevalence, and death became negative after 2000. Tuberculosis in HIV-negative individuals disproportionately occurs in men and boys (versus women and girls); 64·0% of cases (63·6 to 64·3) and 64·7% of deaths (60·8 to 70·3). Globally, malaria cases and deaths grew rapidly from 1990 reaching a peak of 232 million cases (143 million to 387 million) in 2003 and 1·2 million deaths (1·1 million to 1·4 million) in 2004. Since 2004, child deaths from malaria in sub-Saharan Africa have decreased by 31·5% (15·7 to 44·1). Outside of Africa, malaria mortality has been steadily decreasing since 1990. INTERPRETATION: Our estimates of the number of people living with HIV are 18·7% smaller than UNAIDS's estimates in 2012. The number of people living with malaria is larger than estimated by WHO. The number of people living with HIV, tuberculosis, or malaria have all decreased since 2000. At the global level, upward trends for malaria and HIV deaths have been reversed and declines in tuberculosis deaths have accelerated. 101 countries (74 of which are developing) still have increasing HIV incidence. Substantial progress since the Millennium Declaration is an encouraging sign of the effect of global action. FUNDING: Bill & Melinda Gates Foundation.

Global, regional, and national levels and causes of maternal mortality during 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013.

BACKGROUND: The fifth Millennium Development Goal (MDG 5) established the goal of a 75% reduction in the maternal mortality ratio (MMR; number of maternal deaths per 100,000 livebirths) between 1990 and 2015. We aimed to measure levels and track trends in maternal mortality, the key causes contributing to maternal death, and timing of maternal death with respect to delivery. METHODS: We used robust statistical methods including the Cause of Death Ensemble model (CODEm) to analyse a database of data for 7065 site-years and estimate the number of maternal deaths from all causes in 188 countries between 1990 and 2013. We estimated the number of pregnancy-related deaths caused by HIV on the basis of a systematic review of the relative risk of dying during pregnancy for HIV-positive women compared with HIV-negative women. We also estimated the fraction of these deaths aggravated by pregnancy on the basis of a systematic review. To estimate the numbers of maternal deaths due to nine different causes, we identified 61 sources from a systematic review and 943 site-years of vital registration data. We also did a systematic review of reports about the timing of maternal death, identifying 142 sources to use in our analysis. We developed estimates for each country for 1990-2013 using Bayesian meta-regression. We estimated 95% uncertainty intervals (UIs) for all values. FINDINGS: 292,982 (95% UI 261,017-327,792) maternal deaths occurred in 2013, compared with 376,034 (343,483-407,574) in 1990. The global annual rate of change in the MMR was -0·3% (-1·1 to 0·6) from 1990 to 2003, and -2·7% (-3·9 to -1·5) from 2003 to 2013, with evidence of continued acceleration. MMRs reduced consistently in south, east, and southeast Asia between 1990 and 2013, but maternal deaths increased in much of sub-Saharan Africa during the 1990s. 2070 (1290-2866) maternal deaths were related to HIV in 2013, 0·4% (0·2-0·6) of the global total. MMR was highest in the oldest age groups in both 1990 and 2013. In 2013, most deaths occurred intrapartum or postpartum. Causes varied by region and between 1990 and 2013. We recorded substantial variation in the MMR by country in 2013, from 956·8 (685·1-1262·8) in South Sudan to 2·4 (1·6-3·6) in Iceland. INTERPRETATION: Global rates of change suggest that only 16 countries will achieve the MDG 5 target by 2015. Accelerated reductions since the Millennium Declaration in 2000 coincide with increased development assistance for maternal, newborn, and child health. Setting of targets and associated interventions for after 2015 will need careful consideration of regions that are making slow progress, such as west and central Africa. FUNDING: Bill & Melinda Gates Foundation.